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We study 23 to 30 nm long suspended single-wall carbon nanotube quantum dots and observe 
both their stretching and bending vibrational modes. We use low-temperature DC electron transport 
to excite and measure the tubes' bending mode by making use of a positive feedback mechanism 
between their vibrations and the tunneling electrons. In these nano-electro-mechanical- systems 
(NEMS), we measure fundamental bending frequencies fbend ~ 75 - 280 GHz, and extract quality 
factors Q ~ 10 6 . The NEMS' frequencies can be tuned by a factor of two with tension induced 
by mechanical breakj unctions actuated by an electrostatic force, or tension from bent suspended 
electrodes. 



Single-wall carbon nanotubes (SWCNTs) offer the 
prospect of making ultra-high frequency, /, and qual- 
ity factor, Q, nano-electromechanical systems (NEMS) 
due to their small mass, large bulk modulus and tunable 
electro-mechanical properties [TJ |2] . There is a consider- 
able ongoing effort to understand the physics of electron- 
vibron (e — v) coupling for the various vibronic modes 
of carbon nanotubes [3HT3]. An understanding of how 
the strength of these couplings can be controlled with 
mechanical strain [7] could lead to the development of 
NEMS whose vibrational quantum state could be ma- 
nipulated to create mechanical qubits [14] . The SWCNT 
bending mode is of particular interest for potential appli- 
cations, since its frequency, fbend, can be tuned by adjust- 
ing the length and tension of the nanotube, and is very 
sensitive to the tube's environment (adsorbed particles, 
forces). In high-Q devices, the bending mode vibrations 
can couple strongly to electron transport [lOj [TTJ [13J [15] , 
and be used to develop high-resolution mass sensors [TSl - 
HB1 and high-/ oscillators [T9j[20]. Developing SWCNT- 
NEMS whose fundamental mode fbend is beyond the few- 
GHz range remains a major challenge both in terms of 
fabrication, as devices must be a few 10s of nm long, and 
detection, since measuring frequencies in the 10s or 100s 
of GHz can be extremely challenging. 

We report the observation of strain-tunable bending 
mode resonances up to ~ 280 GHz, and Q-factors of the 
order of 10 6 in 23 to 30 nm long suspended SWCNT 
QDs. The e — v coupling for a bending vibron, \bend, is 
much smaller than the ones for the stretching, \ s tretch, 
or breathing, Xbreath, vibrons [7 . However, the effec- 
tive coupling between the tunneling current and bending 
mode is enhanced in high-Q devices, where a large popu- 
lation of bending vibrons can be created due to a positive 
feedback between the electron flow and the vibrations 
[lQj HU [131 [15]. We use DC electron transport in the 
quantum dot (QD) regime, T = 4.2 K unless specified, 
to excite both bending and stretching vibrons, and mea- 
sure their frequencies through their effects on conductiv- 
ity. We demonstrate that fbend can be tuned via strain 
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FIG. 1: (Color online.) Suspended SWCNT quantum dots, 
(a) Top-view SEM image showing an electromigrated gold 
breakj unction across which a 22 nm long suspended SWCNT 
(Device A) is clearly visible (inset), (b) Top view SEM image 
showing a 78 nm long suspended SWCNT. A gold meniscus 
from the left contact extends over the tube and shortens the 
length of the suspended quantum dot (Device C). (c) and 
(d) Two-dimensional I — Vb — Vg Coulomb blockade data 
for Devices A and C respectively. The positive (negative) 
numbers in the Coulomb diamonds refer to the number of 
electrons (holes) in the QD ground state. 



applied to the tube, and extract e — v coupling factors 
^stretch ~ 1, ^bend ~ 10 -3 . These results demonstrate 
that SWCNTs can be used as extremely high frequency 
NEMS, up to a few 100s of GHz, and that their tunabil- 
ity makes them good candidates to study the quantum 
mechanics of macroscopic systems. 

We fabricated several ultra-short, 10 to 80 nm, sus- 
pended SWCNT devices, see Fig. la-b and the Sup- 
porting Information (SI). We observed stretching vibron 
modes in the five devices that we studied in detail, and 
bending mode resonances in three of these devices (De- 
vices A, B, C). We focus our discussion on the later three 
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devices, while data for the other two devices, « 10 nm 
long, are shown in section 1 of the SI. In Fig. la and its 
inset, we see scanning electron microscope (SEM) im- 
ages of the central portion of bow-tie suspended gold 
electrodes (breakj unction) contacting a 22 ± 4 nm long 
section of SWCNT (Device A). Figure lb, shows Device 
C, where a 78 nm long section of suspended SWCNT is 
contacted by rectangular gold electrodes. The complete 
sample fabrication procedure for these devices is detailed 
elsewhere [21], and briefly summarized below. 

To fabricate bow-tie breakj unctions on top of SWC- 
NTs (Devices A and B in main text, and Devices D and 
E in SI), we start with heavily-doped Si wafers with a 
300 nm-thick Si02 film on their top side. The Si sub- 
strate acts as a back-gate electrode. We grow SWCNTs 
by chemical vapor deposition [22] , and use AFM to mea- 
sure their diameter to be d = 2 ± 0.6 nm. We use e-beam 
lithography to define 40- nm thick gold bow-tie junctions 
(no adhesion layer) on top of selected SWCNTs. The 
bow-ties are approximately one micron long and 350 nm 
wide at their center, Fig. la and Figs. S2-S3, and con- 
nected via 3 /im-wide gold wires. We use a wet buffered 
oxide etch (BOE) to remove Si02 from under the gold 
bow-tie junctions to suspend them. The final fabrica- 
tion step is to create a gap in the suspended gold bridges 
by electromigration [2TJ [23H26] to uncover nm-sized sec- 
tions of the SWCNTs. To do so, we ramp up a bias 
voltage, Vb, across the gold bridges and electromigrate 
our devices while they are immersed in liquid Helium or 
in high-vacuum (< 10" 6 Torr) at T « 4.2K. We care- 
fully control the rate of electromigration with a custom 
feedback software to adjust the size of the gap and avoid 
damaging the SWCNTs [21]. The inset of Fig. la shows 
an enlargement of the lower portion of the breakj unction 
of Device A, where the short section of SWCNT across 
the gap is visible. This position matches the location of 
the SWCNT before gold deposition. To fabricate Device 
C, shown in Fig. lb, we directly exposed two rectangular 
gold electrodes on top of the tube using e-beam lithog- 
raphy, and then suspended the tube with a BOE etch. 
Figure lb shows contrast variations along the tube which 
suggests that the gold meniscus from the left contact ex- 
tends over much of the suspended tube, shortening the 
freely suspended section of Device C. 

Device A is a semiconducting tube, and enters the 
QD regime at low-temperature due to Schottky barri- 
ers at the Au/SWCNT interfaces. Figure lc shows the 
low-temperature current, /, versus Vb, and gate volt- 
age, Vg, data for Device A. We observe a clear single- 
electron transistor effect with on and off states, defin- 
ing a Coulomb diamond for each charge ground state of 
the QD [27 . The diamonds are labeled with positive 
(negative) integers, TV, corresponding to the number of 
electrons (holes) on the dot. The wide blockaded region 
around Vq = corresponds to the charge neutrality point 
where no carrier is present on the tube [28] . The location 
of this depletion region near zero gate voltage indicates 
that there is little chemical doping, and that most of the 



contamination adsorbed on the tube during fabrication 
was ashen during the electromigration procedure. From 
the width of the diamonds for N ^ 0, we extract a gate 
capacitance Cq = c/AVq = 0.31 aF [27 . The capaci- 
tance per unit length of the tube can be estimated using 
a wire over a plane model, C/Lq = 27re/(cosh _1 (/i/r)), 
where h is the wire to plane distance and r the wire ra- 
dius. We model Cq as two such capacitors in series (SI 
section 3), respectively with vacuum and Si02 dielectrics. 
For Device A, t ox ide = 100 nm, measured by ellipsome- 
try, and t vac — 35 ± 5 nm as measured by AFM. We 
extract Lq = 27 ± 4 nm. This is in agreement with the 
length L = 22 ± 4 nm measured by SEM in Fig. la, and 
confirms that the capacitance model gives a reasonable 
estimate of the length. Since Lq ~ a single QD occu- 
pies the full length of the exposed SWCNT. Device B was 
fabricated following the same procedure as Device A, but 
its SWCNT is metallic and enters the QD regime at low 
temperature due to imperfect contacts [28]. Its Coulomb 
blockade data, Fig. 2b and SI section 4, show that a single 
QD, whose length is Lq = 30 ± 4 nm, occupies the sus- 
pended tube. The I — Vb — Vq data for Device C in Fig. 
Id, taken at T = 2.2 K, are also consistent with a single 
QD, showing only one set of positive and negative slope 
tunneling thresholds. From the Coulomb diamonds, we 
extract Lq = 23 ± 4 nm for Device C, indicating that 
the QD is shorter than the distance between the gold 
contacts visible in Fig. lb. 

Figure 2 shows dl/dVB — Vb~ Vg electron spectroscopy 
data for Devices A, B and C, allowing us to measure the 
energy of their excited states appearing as resonances 
running parallel to the Coulomb threshold peaks [27] . 
We extract the energy of these excitations as E — eVs 
where the tunneling resonance intersects the Coulomb 
diamonds as indicated by the arrows in Fig. 2a-c. The 
lowest excitation energies are 9, 10 and 7 meV for De- 
vices A, B, and C respectively. These energies are more 
than an order of magnitude too small to correspond to 
electronic excitations of 23-30 nm long QDs [28 . Rather, 
they match the energy scale of the longitudinal stretch- 
ing vibrons [9j [29], AE S = hv s Aq, where Aq = tt/L s , 
L s is the length of the oscillator, and v s = y/Y/p is the 
average group velocity [30]. Our devices have tube di- 
ameters d ~ 2 nm, mass density p = 870 kg/m 3 , and 
Young's modulus Y « 1 TPa, giving v s = 3.4 x 10 4 m/s. 
The extracted length of the longitudinal vibrons in the 
3 devices are respectively L s = 8, 7 and 10 nm. These 
lengths are of the same scale, but shorter than the size of 
the QDs, similarly to what has been reported previously 
in suspended SWCNT devices where the discrepancy was 
attributed to localized vibrons [9j [I2j [31] . 

The excitation energies of the first vibronic sideband, 
extracted at each Coulomb diamond (charge state) be- 
tween Vq = -2 and 7 V for Device B show no signifi- 
cant change (Fig. S4). Vq creates an electrostatic force 
which bends the bow-tie suspended gold contacts and 
can strain the suspended tube by up to a few percent 
(SI section 2). The absence of strain dependence for the 



3 



(a) 40 

20 
£ 

a? 

-20 
-40 




-0.5 0.5 

dl/dV(juS) 



-0.50 -0.25 

V G (V) 



0.0 0.8 

dl/dV (juS) 



3.75 4.00 



(c) 



40 



20 



5 




-20 



-40 



-4.8 -4.4 -4.0 -3.6 -3.2 -2.8 

20 

dl/dV(juS) 



FIG. 2: (Color online.) Stretching vibrons. (a), (b), (c) 
dl/dVB —Vb—Vg data for Devices A, B and C respectively, 
before current annealing. The arrows identify the position of 
some of the vibron sidebands. The red arrows point to side- 
bands which are not always parallel to the Coulomb thresh- 
olds and suggest the emergence of a bending mode resonance. 
In panel (c), the black dashed lines highlight well resolved 
bending mode resonances. 



excitations is consistent with stretching vibrons, as their 
energy is unaffected by strains up to several percent [32] , 
as opposed to bending mode vibrations whose energies 
are strongly strain dependent [33]. Bending mode vibrons 
couple quadratically to electrons and have a much smaller 
e — v coupling than stretching vibrons [7] . For 23-30 nm 
long SWCNTs, bending vibrons have an energy scale of ~ 
0.1 meV which is below our temperature limited energy 
resolution. The radial breathing mode of our SWCNTs 
has a comparable energy scale to the one for stretching 
vibrons, AEb = 28 meV/d where d is the nanotube's di- 
ameter. However, the breathing mode is an optical mode 
and its coupling to electrons is expected to be an order 
of magnitude smaller than for the stretching mode in our 
devices [7]. The breathing mode has only been observed 
in transport experiments when the electrons tunnel in 
the SWCNTs with their momentum perpendicular to the 



tube's axis, as in STM experiments [3j|4]. 

Therefore, we expect that the stretching vibrons cou- 
ple the most strongly to electron transport in our de- 
vices, and we ascribe the excited state sidebands in Fig. 
2 to stretching vibrons. For an electron tunneling on 
a QD while simultaneously creating stretching vibrons, 
the tunneling probability is proportional to the square 
of the overlap of the QD vibronic wavefunction before 
(ground state) and after (excited state) tunneling, and 
dl/dV oc e~ 9 g n /n\, where g — X 2 stretch and n is the num- 
ber of vibrons generated during the tunneling process 
(n = 1 is the first sideband, and so on) [9] [34]. Hence, the 
observation of multiple stretching sidebands with simi- 
lar amplitudes, as in Fig. 2a-b, is only possible for large 
e - v coupling, when \ str etch ~ 1. For \ stre tch > 1, the 
ground state to ground state transition, n = 0, is ex- 
ponentially suppressed at low Vb where the energies of 
the electrons are not sufficient to create vibrons (Franck- 
Condon blockade) [9] [12] [29] EJ • We observe this effect 
in Fig. 2c for Device C, where dl/dV along the Coulomb 
threshold is heavily suppressed near Vb = 0. Close in- 
spection of the positive Coulomb threshold line for the 
N = charge state, shows that its intensity steps up 
around Vb — +7 mV indicating that the n — 1 sideband 
of the stretching mode is present, but suppressed by the 
FC blockade, while the sideband for n = 2 around Vb = 
+15 mV is much stronger. Further evidence of strong 
e — v stretc h coupling, and high Q-factor, is given by the 
nearly horizontal co-tunneling resonance [35] which orig- 
inates where the stretching sideband meets the Coulomb 
threshold at Vb = +15 mV. We conclude that X s tretch is 
larger than 1 in Device C, and comparable to 1 in A and 
B. A small difference in the tubes' diameters may explain 
this variation since X s tretch ^ 1/V^ [7]. 

Figure 2c shows additional features in dl/dV (dashed 
black lines) which are not parallel to the Coulomb thresh- 
old signaling the bending mode resonance of the NEMS. 
The coupling of electrons to bending vibrons, Xbend ~ 
40xL[/im]/(7T(i[nm]) 3 , for the fundamental flexural mode 
is small [7]. Thus, we would not expect to observe ex- 
cited state sidebands in dl/dV when cVb = hujbend as f° r 
the stretching mode, even if our temperature were below 
the bending vibron energy. However, feedback mecha- 
nisms can enhance the effective coupling between bend- 
ing vibrations and tunneling electrons in QDs. Bending 
vibrations produce a time dependent gate capacitance, 
or equivalent ly an effectively oscillating AC gate voltage. 
In the QD regime, this AC gate voltage can create AC 
modulated tunneling rates for the electrons as it moves 
the energy of the dot across one of the Coulomb block- 
ade tunneling threshold. In turn, the electrons tunneling 
on/off the QD create an AC electrostatic force on the 
tube which can further drive the bending mode. When 
the current flow through the QD is large enough so that 
the tunneling rate of electrons coupling to the bending 
vibrons, T, matches or exceeds fbend, and the Q-factor is 
sufficiently high for the vibrons to be long-lived compared 
to the tunneling rate, this can lead to a positive feedback 
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between the bending mode and tunneling electrons which 
spontaneously drives large-amplitude bending vibrations 

From Fig. 2c at V G =-3.9 V and V B = -15.2 mV, 
where the resonance is sharp and extends slightly into 
the normally blockaded region, we extract / = 45 ± 10 
nA and fbend ~ I/ e ~ 280 GHz. We can compare this 
frequency with the expected frequency for a tube under 
zero tension, T = 0, 



(a) 



jw ) r=o = 22.38r/(47ri: 2 ) 1 



(1) 



where r is the tube radius. Using Eq. 1 for Device C with 
L — Lq — 23 nm, we find fbend,T=o = 113 GHz. The 
discrepancy between the measured and expected fbend in- 
dicates that the tube is under considerable tension. The 
predicted dependence of fbend for a SWCNT under high 
tension is [33] 



fbend — 



2tt 




(2) 



where T is the tension. The tension in Device C was 
built-in when the device was suspended (BOE etch), and 
the electrodes contacting the SWCNT bent toward the 
substrate, Fig. S5. Using tilted SEM and AFM to im- 
age the suspended gold electrodes contacting the tube, 
we confirm that the tube is under a considerable tensile 
strain, and measure a strain of 9.8% (SI section 5). This 
strain level is realistic in SWCNTs which can sustain a 
strain up to at least 13.7 % without any plastic deforma- 
tion [36]. Using the measured strain and Eq. 2, the rela- 
tionship fbend ~ I/e predicts I = 36 ± 1 nA (225 GHz), 
where the uncertainty comes from the uncertainty on the 
tube's diameter. This expected fbend agrees with the 
measured current at the resonance, I = 45 ± 10 nA (281 
± 63 GHz), where the uncertainty is due to the sharpness 
of the current step at the resonance. This agreement, to- 
gether with the detailed measurements on Device A dis- 
cussed below, confirms the relationship fbend ~ I/e at 
the bending resonance. We only observed this approx- 
imate equality between fbend an d the current when the 
resonance extends inside the normally Coulomb block- 
aded region. We extract the e — Vbend coupling strength 
in Device C from [37 X be nd = (F 2 ) / '(hm e f f(2n fbend) 3 ) ~ 
10 -3 , where F is the electrostatic force on the tube, and 
m e ff =0.725m tu be is the effective mass of the oscillator. 

We note that the bending resonance in Fig. 2c emerges 
from the second stretching mode sideband at positive Vb 
and from the first stretching sideband at negative Vb- 
This can be explained by the sudden increase of I at the 
vibron sidebands which can trigger the onset of the pos- 
itive feedback necessary for observing the flexural mode. 
The presence of the stetching vibrons also leads to an 
increase in the population of bending vibrons because 
the main decay mode of a stretching vibron is into two 
bending vibrons [38] . 

We do not observe a sharp bending mode resonance 
for Devices A and B in Fig. 2, but we note broad and 
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FIG. 3: (Color online.) Bending resonances, (a) dl / dVu — 
Vb — Vg data for Device A after current annealing. A strong 
non-linearity in dl/dV appears at high Vb (dashed black 
lines) and signals the excitation of the bending mode. Note 
that these new features are not parallel to the Coulomb tun- 
neling thresholds and even extend into the normally blockaded 
region, (b) and (c) One-dimensional I — Vb and dl/dVB — Vb 
data cuts respectively, extracted from the red (thin) and black 
(thick) vertical cuts in (a). 



strong sidebands at high V B (red arrows) which are not 
parallel to the Coulomb diamonds. These features sug- 
gest that the bending mode non-linearity is washed out 
due to a low Q-factor, which would explain the very grad- 
ual rise of / in this region. The current inside the broad 
resonances at the top left and top right of Figs. 2a and 
2b respectively are approximately 15 and 10 nA, which 
correspond to fbend ~ 100 GHz in agreement with the 
expected fbend,T=o fc> r Devices A and B. We used in situ 
high-current annealing to attempt to remove any remain- 
ing surface contamination from Devices A and B which 
may have been limiting their Q. Figure 3a shows the 
same dl/dV — V b — Vq region as Fig. 2a for Device A af- 
ter an annealing step where we flowed 7^4 jiA at V B = 
0.33 V for 10 minutes. A very strong bending mode reso- 
nance, marked by the black dashed line, is now visible at 
positive Vb where previously in Fig. 2a the broad dl/dV 
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feature was present, and an additional resonance is visi- 
ble at negative Vb- For Device B, high current annealing 
at / ~ 4.2 /iA and Vb — 0.66 V for 10 minutes did not 
substantially modify the dl/dV data. The differing ef- 
fectiveness of the annealing procedure on samples A and 
B may be due to the higher resistance of Device B. Since 
the two devices have almost equal lengths, and the tube 
in Device B is metallic, the higher resistance in B is likely 
a result of a higher contact resistance. This would mean 
that a larger portion of the Joule heating power is dissi- 
pated in the contacts, thus the temperature rise of Device 
B during annealing is smaller than in A and the annealing 
less effective in removing contamination. 

There are additional features allowing us to distinguish 
the stretching mode sidebands from the bending mode 
dl/dV resonances in Fig. 3. In panel (a), as was observed 
in Device C, the flexural resonance for Device A extends 
inside what would normally be the Coulomb blockaded 
region. The observed shape of the bending resonances 
in Figs. 3a and 2c are consistent with both theoretical 
calculations [15] and previous experimental observations 
in devices with lower fbend [101 E]. Panels 3b-c show 
respectively 1-dimensional I — Vb and dl/dV — Vb cuts 
along the red (thin) and black (thick) vertical lines in 
panel (a). The width of the bending mode dl/dV res- 
onance is not temperature limited, and can be much 
sharper than both the charge tunneling thresholds and 
the stretching mode side bands. Similarly, the amplitude 
of the change in / due to the bending resonance can be 
much bigger than the one due to the longitudinal mode, 
as seen in the red data in Fig. 3b-c. In Fig. 3a, the bend- 
ing mode resonances show up as positive dl/dV (PDC) 
when they have negative slopes for Vb > 0, or positive 
slopes for Vb < 0. Conversely, dl/dV is negative (NDC) 
when the resonances have positive slopes for Vb > or 
negative slopes for Vb < 0. The same alternating pat- 
tern of (dl/dV) bend between PDC and NDC is observed 
in Fig. 2c for Device C, and in Fig. 4a for Device A over 
a broader range of Vg- While the details of this pattern 
are not currently understood, a possible explanation may 
relate to the QD energy level alignment with the Fermi 
level of the electrodes at the Coulomb thresholds that are 
closest to the bending resonances. A PDC bend mode is 
observed when the resonance is near the Coulomb thresh- 
olds which add an electron to the QD, while a NDC is 
seen near the transitions which remove an electron from 
the dot. Thus, the bending vibrons modulate the tunnel- 
ing rate on (off) the QD at PDC (NDC). Since vibrons 
are created by electrons tunneling on the QD, at PDC 
peaks the vibrons are in phase with the incoming elec- 
trons and can help scatter them across the dot. On the 
other hand, at the NDC peaks the vibrons modulate the 
rate of electrons tunneling off the dot which may cause 
an out-of-phase driving of the vibrons and lead to more 
backscattering of the electrons. 

The 1-dimensional vertical red (thin) cut in Fig. 3a in- 
tersects the bending resonance at negative Vb in a region 
which would normally be inside the Coulomb diamond, 



and where the relation fbend ~ T « I/e holds. Using a 
Lorentzian fit, we extract I = 12 ± 6 nA at the dl/dV 
resonance for the red (thin) data in Fig. 3b-c located 
at V B = -38.3 mV and V G = -0.55 V, giving f bend w 
75 ± 30 GHz. Using Eq. 1, we calculate for Device A, 
fbend,T=o = 84 GHz. This close agreement suggests that 
Device A is not under a significant tension at Vq = -0.55 
V. The expected strain from the electrostatic bending of 
the suspended bow-tie cantilevers due to Vq = -0.55 V 
is negligible, and SEM and AFM imaging of Device A 
show no apparent static (built-in) bending of the elec- 
trodes which would strain the SWCNT (SI section 2). 
To calculate the Q-factor of Devices A and C based on 
their transport characteristics, we follow Lassagne et al. 
[TTj [37] who derived for similar samples, 



1 _ 27T fbend ( 2CqVq_ 



G 



(3) 



where C' G is the derivative of Cq with respect to the verti- 
cal displacement of the tube, k is the effective spring con- 
stant of the tube, is the total capacitance of the QD, 
and G the conductance (SI section 6). Using data from 
Figs. 3a and 2c for the bending resonances located inside 
the normally Coulomb blockaded region, we calculate for 
Device A, at Vq =-0.55 V, and Device C, at Vq = -3.9 
V, respectively Q = 2.9 xlO 6 and 1.3 x 10 6 . These high 
Qbend values are expected [38] since Qbend 1/L 2 , and 
to the best of our knowledge, our devices are the short- 
est SWCNT-NEMS reported thus far. The calculated 
Q-factors are also consistent with previous work where 
bending mode resonances were observed in DC transport 
only when Q > 10 4 [101CE3]. 

Figure 4 shows data for Device A over a broad range 
of Vg- The black markers indicate the locations of the 
bending dl/dV resonances where we extract the / val- 
ues shown in Fig. 4b. The positive bias dl/dV bend- 
ing resonances are left unmarked to show the alternat- 
ing PDC/NDC pattern at each Coulomb diamond. As 
Vg increases, the electrostatic force between the gate 
electrode and the suspended bow-tie gold electrodes in- 
creases oc Vq. Using the geometry of the two suspended 
cantilever contacts making up the bow-tie junction we 
can calculate the vertical displacement of the suspended 
source and drain electrodes, Ay, and the longitudinal 
strain AL/L, induced by the moving contacts (SI sec- 
tion 2). This strain applies a tension to the tube which 
increases rapidly with Vg, and modifies fbend as given by 
Eq. 2. 

Fig. 4b shows I extracted from the negative Vb reso- 
nance in panel (a) versus Vg- The current oscillates with 
the same period as the Coulomb diamonds. The data in 
Fig. 4a has coarser steps in Vg than in Figs. 2 and 3, 
which limits the resolution of the / oscillations in Fig 4b. 
The effect of the width (resolution) of the Coulomb dia- 
monds on the depth of the current minima is made clear 
when comparing the data from the far left and far right 
of Fig. 4a-b. While the high positive gate diamonds are 
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FIG. 4: (Color online.) Tension-tuning of fbend- (a) 
dl/dVB — Vb — Vg data for Device A, after current anneal- 
ing, over a broad range of Vg- The gate voltage bends the 
breakj unction contacts and strains the suspended SWCNT, 
which tunes fbend- The black markers indicate the locations 
where we extracted the bending mode resonance I values 
shown in (b) . (b) / measured at the bending resonances with 
Vb < versus Vg- The red (thin) and blue (thick) lines 
show the low and high-tension regime theoretical calculation 
for / = e fbend valid at the minima of well-resolved Coulomb 
valleys. Due to the coarseness of the Vg data steps, the cal- 
culated current is only expected to agree quantitatively with 
the data at the best resolved (widest) Coulomb valleys located 
around Vg = and -4.2 V. 

narrow and their / oscillations poorly resolved, the neg- 
ative Vg data show wider diamonds and deeper current 
minima. The Coulomb valleys around Vg = V and -4 
V are wide and show the best denned current minima, 
where we expect / « e fbend- The red (thin) line and 
blue (thick) line in Fig. 4b are respectively the calculated 
low-tension regime (essentially Eq. 1) and high-tension 
regime (Eq. 2) current expected when the bending res- 
onance is inside the Coulomb blockaded region. These 
calculations do not include any fitting parameters, and 
are based on the measured dimensions of Device A (SI 



section 2). The agreement between the data and theoret- 
ical estimate around Vg = V and -4 V, combined with 
a similar agreement for Device C, confirms that we can 
measure fbend with a reasonable accuracy in DC trans- 
port experiments. The frequency around Vg = is « 80 
GHz, and increases up to 150 GHz around Vg = -4.2 V 
where the tube's strain extracted using Eq. 2 is 6.0 %. 

In summary, using DC electron transport in 23-30 nm 
long suspended SWCNT QDs we observed both their 
stretching and bending vibrational modes, and e — v cou- 
plings, Xstretch ~ 1 and X bend -TO -3 . The bending mode 
dl/dV resonance is due to a positive feedback mecha- 
nism between electron tunneling and the bending vibra- 
tions which turns on when the tunneling rate matches 
or exceeds fbend- Using current annealing, we removed 
contamination absorbed on one of the oscillators and 
greatly enhanced the dl/dV resonance associated with 
the bending mode. We measured fundamental bending 
mode frequencies up to ~ 280 GHz, roughly two orders of 
magnitude higher than previously reported fundamental 
frequencies [2j [TUJ [TTJ [19] , and almost an order of mag- 
nitude higher than the highest fbend previously reported 
[20] , The calculated Q-factors of our NEMS are -10 6 . 
The flexural frequency fbend can be tuned by a factor of 2 
with tension from an electrostatically-actuated mechan- 
ical breakj unction, or by a built-in tension due to bent 
suspended contacts. These extremely high fbend and Q- 
factor NEMS are candidates to create extremely sensitive 
force and mass sensors, and to explore the quantum me- 
chanics of macroscopic systems [14] . We plan on probing 
the ultimate limits of SWCNT-NEMS by building even 
shorter suspended QDs, ~ 5-10 nm [2T] . with indepen- 
dently tunable charge states and mechanical strain using 
a gated-mechanical breakj unction approach [25] . 
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